Muons from Neutralino Annihilations in the Sun: Flipped SU(5) 
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We consider two classes of supersymmetric flipped SU(5) models with gravity mediated supersymmetry 
breaking such that the thermal neutralino relic abundance provides the observed dark matter density in the 
universe. We estimate the muon flux induced by neutrinos that arise from neutralino annihilations in the Sun and 
discuss prospects for detecting this flux in the IceCube/Deep Core experiment. We also provide comparisons 
with the corresponding fluxes in the constrained minimal supersymmetric standard model and non-universal 
Higgs models. Regions in the parameter space that can be explored by the IceCube/DeepCore experiment are 
identified. 

PACS numbers: 95.35.+d, 98.62.Gq, 98.70.Vc. 95.85.Ry 



I. INTRODUCTION 



It is now widely accepted that approximately 23 % of the Universe's energy density consists of non-baryonic cold dark matter 
[1]. A large number of experiments consisting of direct, indirect and accelerator searches are currently underway all hoping 
to discover the underlying, presumably massive (~ GeV -TeV), weakly interacting dark matter particle (WIMP). The lightest 
i-C ' neutralino in supersymmetric models with conserved matter parity is a particularly attractive cold dark matter candidate and 
^p-i has attracted a great deal of attention. The direct detection searches have already yielded important constraints on the spin 
independent neutralino-nucleon cross sections in the constrained minimal supersymmetric standard model (CMSSM) and some 
D [ related models (see [2], and references therein). 

^ ' Indirect WIMP searches rely on the capture and subsequent annihilation, say in the Sun's center, of relic dark matter particles. 
[ The neutralinos, in particular, can annihilate into the known SM particles, for example, XX t'^t~ . The tau particles in 
fSJ , turn produce energetic muon neutrinos which interact with the polar ice to produce muons which can be identified by the km^ 
I> ' IceCube/Deep Core detector [3]. 

. Neutralinos in the galactic halo passing through a massive body like the Sun can get captured if they scatter off the nuclei 
' with velocities smaller than the escape velocity. In the core of the Sun, where they eventually accumulate, these neutralinos can 
annihilate into known SM particles, for e.g., xX ~^ t~^t~ . These particles decay (e.g. t — > ^^T-i/^/i) and produce energetic 
muon neutrinos which can then be detected at IceCube after they interact with the polar ice and produce muons (for e.g. via 
processes like v^^ + N ^ /i^+X,N being the nucleon and X some hadronic system). We investigate the possibility of detecting 
these energetic neutrinos by estimating the flux of muons that they induce. In addition to the well studied CMSSM, we explore 
other well motivated models, namely, flipped SU(5), non-universal Higgs models (NUHM2), and flipped SU(5) with universal 
Higgs masses at Mgut- The prospects for detecting this neutrino induced muon flux by the IceCube/DeepCore experiment is 
discussed. 

In this paper we are mainly interested in studying the implications of supersymmetric flipped SU(5) models for indirect dark 
5J] matter WIMP searches, with the lightest neutralino being the dark matter candidate. Flipped SU(5) has several distinct features 
^ which are not easily replicated in other GUTs such as SU(5) and SO(IO). For instance, the well-known doublet-triplet splitting 
problem is easily solved in flipped SU(5) [4]. Primordial inflation with predictions for the cosmological parameters in good 
agreement with the 7 year WMAP data are readily obtained , which in turn, lead to testable predictions for proton decay [5]. 

Our paper is organized as follows. In section II, following [6], we briefly describe the two flipped SU(5) models under 
discussion. Consistent with the underlying gauge group (SU(5) x U(l)), both classes of models work with non-universal gaugino 
masses. Their difference stems from the non-universal soft scalar Higgs masses employed in one of the models. In section 
III we review the calculations of the conversion factors relating the muon flux and spin dependent (SD) cross section in the 
IceCube/Deep Core experiment. Section IV contains a description of the experimental constraints and the scanning procedure 
employed to generate the benchmark points. Our predictions for the muon flux and SD cross section are presented in section V, 
and the conclusions are summarized in section VI. 
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II. THEORETICAL FRAMEWORK OF THE MODELS 

We are interested in estimating the neutrino induced muon fluxes from neutralino annihilations in the Sun, with flipped 
SU(5) (FSU(5)) boundary conditions imposed on the soft supersymmetry breaking (SSB) parameters at Mqut- More generally, 
we compare four distinct models, namely CMSSM, NUHM2, FSU(5) and FSU(5) with Universal SSB Higgs mass boundary 
condition (FSU(5)-UH). We will briefly describe each model below. The CMSSM [7] has the following parameters at the 
Mgut- 

Too, »7ii/2,Ao, tan/3, sign(/x). (1) 

Here toq is the soft supersymmetry breaking (SSB) scalar masses, toi/2 is the SSB gaugino mass, is the universal SSB 
trilinear scalar interaction (with the corresponding Yukawa coupling factored out), tan/3 is the ratio of the vacuum expectation 
values (VEVs) of the two MSSM Higgs doublets, and the magnitude of fi, but not its sign, is determined by the radiative 
electroweak breaking (REWSB) condition. 

Whereas, universal scalar masses are motivated to suppress unwanted flavor changing neutral currents, the Higgs mass pa- 
rameters can be non-universal. The effects of this non-universality on the parameter space has been studied in models called 
Non-Universal Higgs Models (NUHM) [8]. One of the types of these models called the NUHM2 has two additional parameters 
compared to the CMSSM 

(2) 

where to|^^ and are the SSB the MSSM Higgs mass^ term. The supersymmetric flipped SU(5) (FSU(5)) model [9] is 
based on the maximal subgroup G = SU(5) x U(1)_y of SO(IO), and the sixteen chiral superfields per family of SO(IO) are 
arranged under G as: lOi = (d^, Q, v"^), 5_3 = (w'^, L), I5 = e"^. Here the subscripts refer to the respective charges under 
and we follow the usual notation for the Standard Model (SM) particle content. The MSSM electroweak Higgs doublets 
Hu and Hd belong to 5^ and 5^ of SU(5), respectively. We will assume for simplicity that the soft mass'^ terms, induced at 
7\/gut through gravity mediated supersymmetry breaking [10], are equal in magnitude for the scalar squarks and sleptons of 
the three families. The asymptotic MSSM gaugino masses, on the other hand, can be non-universal. Due to the FSU(5) gauge 
structure, asymptotic SU(3)^ and S\]{2)^ gaugino masses can be different from the U(l)y gaugino mass. Assuming SO(IO) 
normalization for the hypercharge generator in FSU(5) is given hyY = (— ^5/2 + \/2iX)/b, where and X are the 

generators of SU(5) and [11]. We then have the foUowing asymptotic relation between the three MSSM gaugino masses: 

1 24 

Afi = —A/5 + — M', with Ms = Ah - Af3, (3) 
25 25 

where M^, M' , M3, M2 and Mi denote SU(5), U(l)_^-, SU{3)c, SU{2)l and U{1)y gaugino masses respectively. The su- 
persymmetric FSU(5) model thus has two independent parameters {M2 = M3, M') in the gaugino sector. In other words, 
in FSU(5), by assuming gaugino non-universality, we increase by one the number of fundamental parameters compared to the 
CMSSM 

We will also consider both universal {nijj ~m^^) and non-universal (m^ ^rnj^J soft scalar Higgs masses in FSU(5), using 
the notations for FSU(5)-UH and FSU(5) respectively. This would mean up to three additional parameters compared to the 
CMSSM. This latter case, with one additional gaugino mass parameter and two soft scalar mass parameters, provides us with a 
compelling neutralino dark matter candidate for indirect and direct detection [6] in the ongoing and future experiments. 

We use the (fi, m^) parameterization to characterize non-universal soft scalar Higgs masses rather than {Hu,Hd). The funda- 
mental parameters of our FSU(5) model are 

mo, AT', Ar2,tan/3,Ao,At,mA, (4) 

we will assume that /i > 0. Note that /i and ruA are specified at the weak scale, whereas the other parameters are specified at 
A/gut- Although not required, we will assume that the gauge coupling unification condition 93 = gi = 92 holds at A/qut in 
FSU(5). Such a scenario can arise, for example, from a higher dimensional theory [12] after suitable choice of compactification. 



HI. MUON FLUX, SD CROSS SECTION AND CONVERSION FACTORS 

In this section we will review the calculation of the muon flux and spin dependent cross section and revisit the way the 
conversion factors between the two are calculated. The IceCube collaboration [13] has presented its results as a future bound on 
the muon flux from the Sun. The bound was then converted to a bound on the spin dependent cross section by suitable conversion 
factors calculated in [14] and also discussed in [15]. 
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The flux of neutrino induced muons from neutralino annihilation in the Sun is given by 

Here Ta is the annihilation rate, n is the target number density, and Dq is the distance from the center of the Sun to the detector. 
diV/ /dEi, is the differential energy spectrum of the number of neutrinos from neutralino annihilation with the corresponding 
branching fractions Bf. Ejj^ is the threshold energy of the muon in the detector, A is the muon range, P^u^yi^Efj^^ is 

the survival probability for a muon, dai,{Ei,, E^)/dE^ is the differential neutrino cross-section and P{n, i) is the oscillation 
probability for a neutrino of flavor i to oscillate to flavor /i in the detector. 

The annihilation rate for the weakly interacting massive particles (WIMPs) in the center of the Sun is given by 

Ta = icctanh2(t/r), (6) 

T = {CcCa)~^^^ is a measure of the time in which capture and annihilation equilibrate, Cc is the capture rate, and Ca 
parameterizes the annihilation rate of the WIMPs. For present WIMP annihilation rate, t is the age of the Sun, i.e., f = ~ 
4.5 • 10^ years. The annihilation and capture rate are in equihbrium when t®/r ^ 1, which imphes 

Ta = ^Cc. (7) 

Since the capture rate is proportional to the spin dependent and spin independent cross sections, there would be a direct cor- 
relation between the flux and the SD cross section. The converted bound cannot be trusted for models where the equilibrium 
condition is not satisfied. 

Accurate expressions for the capture rate can be found in [16], while reference [17] gives the approximate expressions. For 
the case of SD cross section, which occurs mainly on hydrogen and the form factor suppression is negligible, the capture rate in 
the Sun can be written as 



r, , 9^ 1, /270kms-i\ / \ /100GeV\ / (Tsd \ , , x 

where ctsd is the neutralino-proton spin dependent cross section, v = 270 Kms~^ is the dark matter velocity dispersion, 
= 0.3 GcV cm~'^ is the local dark matter density, and S{m^/mp) is the kinematical suppression factor defined as 



AV2 



2/3 



with 



{vesc} denotes the mean escape velocity from the Sun. 

Reference [14] calculates accurate conversion factors including neutrino oscillations. Here we take a simple example to see 
how the conversion factors are calculated to get the SD cross section from the muon flux. Ignoring detector thresholds and taking 
the effective range of muons in the detector, the rate of neutrino induced through going muons for the Sun can be approximated 
as [17] 

r, « (1.27X io-^WV-^)pr (y^)'E«'^'E%(^-')/.'(-x)- (11) 

i F 

Since the capture in the Sun is mainly through spin dependent scattering, we can assume asi = {Cc = Cfjj) to get a bound 
on the SD cross section, are the neutrino scattering coefficients = 6.8 and ap = 3.1, and 5; are the muon range coefficients 
with bu = 0.51 and 6p = 0.67. The quantity {Nz^) f_i is the second moment of the neutrino spectrum of type i from final state 
/, scaled by the square of the injection energy Ein of the annihilation products, and is given by 



El J \dEj 



{Nz')f4m^)^-^ ( — ] {E,,E„,)EidE,. (12) 
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The neutrino spectrum from the W~^W and r+r channels can be taken as [18] 

§-Y =^^(l + i?.n) — , (13) 

/ WW 

with E,n{l - /3/2) < < E.nil + 13/2), Tw^^^u = 0.105 , /3 = (1 - ml,/Eiy^\ = 1.01 x IQ-^GcV-i and 
Tp = 3.8 X lO-^GeV-i, and 

dN\ ° ^ 2r^^^ _ ^ 2a;3)(l + (14) 

with < E'l/ < X = E^/Ein and Ft-^^j/i, = 0.18, a^, ~ 5.1 and ap = 9.0. Note that improved functions for the spectra 
can be obtained by using programs like Pythia. In Fig. 1 we show plots of the second moment of these functions. We take Bj = 1 
and only consider contributions of the hard channels W^W~ and r+r" . In our plots we assume that only the W^W^ channel 
contributes for > 80 GcV, and the r+r^ channel for < 80 GcV. The second moments {Nz'^)ww ^nd {N z'^)^.,- for the 
W^W~ and r+r^ channels are obtained by inserting Eqs. (13) and (14) in (12). From Eq. (11), 



C® /TO \ 2 

F^ = (1.27 X 10-23km-V"i)-^ (tg^v) \^-^UN z'')wwA^x) + 'i.Q^N z-')wwA^x)l (15) 



for TO^ > (<)80 GcV for the W+W- (r+r") channel. 



•■X 

We can re-write Eq. (8) as 



and inserting this in Eq. (1 1) we have 



which yields the conversion factor 



CId = /l("^x)<^SD, (16) 
= /i(»7tx)/2(m-x)'^SD, (17) 



^ = ./r\m^)/2-i(TO-,)(W yr cm^). (18) 

The conversion factor and the converted future IceCube/DeepCore bound obtained from it are plotted in Fig. 2. The conversion 
factor we find is similar to the estimates made in [ 14] without including neutrino oscillations. We will see later how the converted 
IceCube/DeepCore future bound for the SD cross section is changed once neutrino oscillations are included. 

Neutrinos created from neutralino annihilations in the Sun's center may experience oscillations and interactions (neutral as 
well as charged) in the Sun. These and other factors, like the Solar composition, can modify the flux observed by detectors 
on Earth. On the other hand, direct detection involves measuring the scattering cross section of neutralinos off of nuclei. The 
conversion of the limit obtained from indirect detection experiments thus involves several uncertainties [14]. They include: 

• Uncertainties in the Solar model. 

• Gravitational effects of planets like Jupiter. 

• Form factors. 

• Variations in local dark matter density and velocity distributions. 

• Neutrino oscillations. 

These uncertainties can affect the estimated muon flux and thereby the deduced cross section. Gravitational affects from Jupiter, 
for example, can reduce the estimated muon flux, whereas neutrino oscillations can enhance it. The form factor suppression is 
negligible for the case of spin dependent interaction since the capture in the Sun mainly occurs through scattering on the hydrogen 
nuclei. All of the above listed affects can have implications for the particular particle physics model being investigated. The 
implications for our models are discussed in section V. 



IV. PHENOMENOLOGICAL CONSTRAINTS AND SCANNING PROCEDURE 



We use the package DarkSUSY-5.0.5 [19] to calculate the flux of neutrino induced muons from the Sun. DarkSusy (DS) uses 
a local dark matter density of = 0.3 GeV/cm"^. From the three methods DS employs to calculate the relic density, we pick 
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the one which includes coannihilations only if the mass difference between the LSP and NLSP is less than 30%. We rescale the 
neutralino density to p = po(ilft.'^/0.025), for flh^ < 0.025. This is done because the neutralino cannot make up all of the dark 
matter in the galaxy halos if Uh^ drops below 0.025. 

Although not required, we assume for simplicity, that the gauge coupling unification condition gi = 92 = .93 holds at Mqut 
for FSU(5). DarkSUSY-5.0.5 uses Isajet 7.78 [20] for renormalization group evolution (RGB) running and the latter employs 
two loop MSSM RGBs and defines Mqc/t to be the scale where gi = 52- A few percent deviation from exact unification 
(51 = 92 = .93) can be attributed to unknown GUT-scale threshold corrections [21]. 

For the random scan we employ the following ranges for our parameters 

< mo < 5 TcV, 
< M' < 1 TcV, 
< M2 < 1 TcV, 
< < 1 TeV, 
< M < 10 TeV, 
tanP = 10, 30, 50, 

^0 = 0, (19) 

we set rrit = 173.1 GeV. 

The random scan is performed over the parameter space of CMSSM {M2=M^=M'=mii2, ■mo='mH^=iTiHj), FSU(5), 
NUHM2 {M2=M'), and FSU(5) with universal soft Higgs masses^ ("'^li^-^'H^) Mqut- Since the neutralino mass is 
sensitive to the gaugino masses, we manipulate the latter to obtain more allowed points in the parameter space. We take piece- 
wise intervals [0,10], [10,100] and [100,1000] for the gaugino masses (in units of GeV). The random points in each of these 
intervals are distributed logarithmically. These points were then combined with a uniform distribution of points on the interval 
[0,1000], with the total number of points around a million, which enables us to obtain a sufficiently dense set of points for our 
plots. This is still not sufficient for FSU(5)-UH, so we perform a Gaussian scan around the allowed points. The code makes a 
Gaussian distribution of points for the scalar and gaugino mass parameters around a point satisfying all the imposed constraints, 
with the variance and mean of the Gaussian distribution being 1/25 and 1 respectively. The random function RNORMX(), 
available in the program library of CBRN, was used to make this Gaussian distribution. 

We apply the experimental constraints on the data sequentially, with all of the collected data points satisfying the requirement 
of radiative electroweak symmetry breaking (RBWSB), and the neutralino in each cases being the LSR On this data, we impose 
the following constraints: 



m~± (chargino mass) > 103.5 GeV [22] 

TOf (stau mass) > 105 GeV [22] 

rUg (gluino mass) > 250 GeV [22] 

TOj (stop mass) > 175 GeV [22] 

TO^ (sbottom mass) > 222 GeV [22] 

m/i (lightest Higgs mass) > 114.4 GeV [23] 

BR{B, fJ-+fJ- ) < 5.8 X 10"^ [24] 

2.85 X 10-'' < BR{b 57) < 4.24 x lO"'' (2cr) [25] 

^cDuh^ = o.iutE;;;^? (5^) [i] 



Note that we do not include the (g — 2)^ constraint for the rest of our discussion. 

V. RESULTS 

We next present the results from the scan over the parameter space listed in Bq. (19). In Fig. 3, we show how the converted 
IceCube/DeepCore future bound is altered with the inclusion of neutrino oscillations. The colored points are consistent with 
RBWSB and satisfy the WMAP relic density bound in the 5(t range, particle mass bounds, and all constraints coming from 
the B-physics. We used different color coding to distinguish different channels for neutralino dark matter. We see from the 
right panel in Fig. 3, and as noted in reference [14], the inclusion of neutrino oscillations has a notable affect on the bound, 
especially for neutralino mass less than the W boson mass. As the W boson decouples (m^ < 80 GeV) and the contribution 
from the tau channel becomes relevant, the bound changes notably. The reason for this is that the tau neutrinos from the decay 
/i^Ppi^T- can oscillate into muon neutrinos, thereby enhancing the muon flux in the detector. This is especially relevant 
for the low neutralino masses we have in FSU(5). In the left panel the future IceCube/DeepCore bound (solid line) indicates that 
the light neutralino {m^ < 70 GeV) parameter space can be tested at IceCube/DeepCore detector, but this same region of the 
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parameter space, it seems, does not yield sufficient muon fluxes, as can be seen in the right panel of Fig. 3. So we see that there 
is some discrepancy. 

Note that, whereas the calculation of the flux is sensitive to the various channels from which the neutrino arises, the cross 
section is not. Fig. 3 shows that as the neutralino mass falls below the W mass, and only soft channels (e.g. bb) are left , the 
muon flux starts decreasing, whereas the SD cross section decreases much less rapidly [26]. This is because the SD cross section 
is not prone to the hardness or softness of the channel, which is not the case for the bound on the SD cross section where the 
sensitivity of the flux is translated to the SD cross section. 

For the calculation of the flux we use the 'default' method in DS which uses the approximate expression for the capture rate 
in the Sun from [17]. It is understood that the dark matter prediction is no longer a natural consequence of supersymmetry [27], 
but requires special relations among the parameters. To have the correct relic dark matter abundance, we require coannihilation, 
resonance or specific Bino-Higgsino mixing solution. On the other hand this then yields some very specific structure for the 
sparticle spectroscopy which can be tested at the LHC. This explains why, in Figs. 4, 5 and 6, we show the various relic channels 
for neutralino dark matter. The colored points are consistent with REWSB and satisfy the WMAP relic density bound in the 5 a 
range, particle mass bounds, and all constraints from B-physics. Figs. 4, 5 and 6 present the muon flux induced by the neutrinos 
originating from annihilating neutralino dark matter in the center of the Sun, for tan/3 = 10, 30, 50. The points shown satisfy 
the WMAP relic density bounds in the 5a range. The calculated muon flux is integrated above a threshold energy of 1 GeV. 
From Figs. 4, 5 and 6, we observe that the IceCube/DeepCore detector can test the following neutralino dark matter scenarios: 
Bino-Higgsino dark matter, light Higgs resonance and finally the "non-identified channel", which is a combination of various 
channels. The points we designate as being non-identified means that the conditions we apply on the neutralino to be from all 
other channels are not satisfied. Thus an observed signal at the IceCube/DeepCore detector can narrow the probable neutralino 
dark matter candidates, and combining this with a signal from the LHC may help identify the nature of dark matter. 

In Fig. 7 we show the results in the fundamental parameter planes. Here Mif2 stands for the GUT scale universal gaugino 
mass in CMSSM and NUHM2, and A/i is the Bino mass for FSU(5)-UH and FSU(5) plots, the expression for which is given in 
Eq. (3). The green points are consistent with REWSB and satisfy the WMAP relic density bound in the 5 a range, particle mass 
bound, and all constraints coming from B-physics. The red points are a subset of the green ones and give muon fluxes which can 
be tested at the IceCube/DeepCore experiment. The IceCube/DeepCore experiment, we see, can test a significant region of the 
flipped SU(5) parameter space. 

In Figures 8, 9, 10 and 11 we show the results in tua vs. rriy^, nif vs. rriy^, vs. and rn^± vs. planes respectively. 
The color coding in these figures is the same as for Fig. 7. We can see that the FSU(5) model gives rise to signals, corresponding 
to a relatively light f and itia, which can be seen by the IceCube/DeepCore experiment. This is not the case for CMSSM and 
NUHM2 models. 

For tan/3 = 10 the lower mass bounds on the lightest neutrahno for the four models (CMSSM, NUHM2, FSU(5)-UH and 
FSU(5)) are 76.7 GeV, 53.1 GeV, 32.2 GeV and 31.6 GeV respectively. This is consistent with a recent study [28] which found 
a lower bound of 28 GeV on the mass of the LSP neutralino. From Fig 5 we see that this neutrahno in the CMSSM case comes 
from the focus point region and when the Bino or Higgsino mixing is large. For NUHM2, it is from the h-resonance channel, 
whereas for FSU(5)-UH it is from the "non-identified" region. We note that the muon flux is highest when the LSP dark matter 
neutralino is mainly Bino-Higgsino like, and this observation is valid for the CMSSM as well as for its extensions. 

Note that the IceCube/DeepCore bound is a conservative one and the muon flux limits can be improved by an order magnitude 
[29]. As noted in [29] their are prospects of detecting the CMSSM focus point (FP) region in the IceCube/DeepCore experiment. 
This can also be seen in Figures 4-6 for the CMSSM. Going from CMSSM to FSU(5) changes this, and in addition to more 
points in the FP region we also have some points from the non-identified and h-resonance regions. The non-universality of the 
Higgs mass parameters opens up the A-funnel region where resonant annihilation occurs through the CP odd Higgs boson A. The 
Bino-Wino coannihilation [30] channel arises from the gaugino non-universality in Eq. (3) and occurs for 2M2 ^ 7\fi at Mqut- 
As we can see from the muon flux plots, this region of the parameter space is not detectable with the current IceCube/DeepCore 
experiment. 

Finally, in TABLE I we present three FSU(5) benchmark points for tan/3 = 10, 30, 50 which yield observable muon fluxes. 
The first point belongs to the stau coannihilation region, the second point is associated with the Bino-Higgsino dark matter with 
light charginos, and the third corresponds to the h-resonance region. 



VI. CONCLUSIONS 

We have considered indirect neutralino dark matter detection in two sets of supersymmetric flipped SU(5) models. These 
two sets of models have non-universal soft gaugino masses at Mqut that are related by the underlying SU{5) x U{l)x gauge 
symmetry. The supersymmetry breaking soft Higgs masses^, associated with i7„ and Hd, are equal at Mqut in one set of 
models (FSU(5)-UH) but not in the other (FSU(5)). We have provided estimates of the flux, from annihilating neutralinos 
in the Sun, of neutrino induced muons, and considered prospects of detecting this flux in the IceCube/DeepCore detector 
Some uncertainties arise in converting the muon flux into spin dependent neutralino-nucleon cross sections that we have briefly 
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discussed. We offer comparisons with previously studied CMSSM and NUHM2 models, and also highlight some benchmark 
models in flipped SU(5) with varying neutralino compositions which can be tested by the IceCube/DeepCore experiment. Our 
results for NUHM2, FSU(5)-UH and FSU(5) show more points above the projected IceCube limit compared to the CMSSM, 
and hence a greater prospect of detection. This is to be expected because the models are less constrained than the CMSSM and 
possess additional free parameters. 
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FIG. 1. Second moment of the neutrino spectrum. The solid lines correspond to neutrinos and the dashed lines are for anti-neutrinos. The 
W^W chaimel is in red and t^t~ channel is in blue. 
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FIG. 2. The left panel shows the conversion factor calculated from approximate expressions for the flux and spin dependent cross section. The 
right panel shows the converted future IceCube/DeepCore muon flux bound using this conversion factor. 




FIG. 3. Comparison of the future IceCube/DeepCore muon flux bound and the converted SD cross section bound. The muon flux is from 
the Sun above 1 GeV threshold is shown for tan/3 = 10. The dashed line in the left panel is the future IceCube/DeepCore bound obtained if 
neutrino oscillations are not included in the flux calculation. The conversion factors used are given in reference [14]. 
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FIG. 4. Flux of fj, + Jl from the Sun above 1 GeV threshold for tan/? = 10. The black line shows the future IceCube/DeepCore bound 
[13]. The colored points are consistent with REWSB and satisfy the WMAP relic density bound in the 5 a range, particle mass bound, and all 
constraints coming from B-physics. The points in different colors correspond to the various solutions of LSP neutralino to be a dark matter 
candidate. 




FIG. 5. Flux of /i + /i from the Sun above 1 GeV threshold for tan/? — 30. We use the same color coding as in Fig. 4. 
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FIG. 6. Flux of fj. + jj. from the Sun above 1 GeV threshold for tan/? — 50. We use the same color coding as in Fig. 4. 
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FIG. 7. Plot in the mi/2 and AIi vs mo plane. We are comparing the allowed parameter spaces for differed models. The green points are 
consistent with REWSB and satisfy the WMAP relic density bound in the 5 a range, particle mass bound, and all constraints coming from 
B-physics. The red points are a subset of the green ones and can generate detectable muon fluxes at the IceCube/DeepCore detector. 
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FIG. 8. Plot in the rriA — rn^ plane. Color coding same as in Fig. 7 
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FIG. 9. Plot in the ruf — plane. Color coding same as in Fig. 7 
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FIG. 10. Plot in the m^^ — plane. Color coding same as in Fig. 7 




FIG. 11. Plot in the m ± — plane. Color coding same as in Fig. 7 
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TABLE I. Sparticle and Higgs masses (in GeV), with rtit — 173.1 GeV. These benchmark points satisfy all the constraints imposed in 
Section V. 



